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Abstract 

The recent results in QCD at low energies are reported. The theoretical 
analysis of hadronic r-decay is performed in complex g^-plane. The terms 
of perturbation theory (PT) up to are accounted, the terms of operator 
product expansion (OPE) - up do dimension 8. The momenta of charmed 
quark vector current polarization operator are analyzed in 3 loops with ac- 
count of dimension 4 and 6 OPE terms. The main conclusion is that in the 
range of low-energy phenomena under consideration PT and OPE are in good 
agreement with experiment starting from = —q^ ~ 1 GeV^. The values of 
as^m^), Usigq)'^ , ((as/7r)G^) and mc(mc) were found. 

1 Introduction 

The first analysis of QCD at low energies in framework of PT and OPE was per- 
formed long ago. The most powerful approach is the method of QCD sum rules 
[1] (for review of old works see [2], more recent results are reviewed in [3]). In the 
earlier analysis it was accepted [1], that QCD coupling constant is small even 
at low energies as{l GeV^) ~ 0.3 (Ag^^, ~ 100 MeV) and only first order as cor- 
rections were accounted. Now it is clear, that at low ~ 1 — 5 GeV^ asiQ"^) is 
about two times larger, then it was taken in [1]. From the other side, new, more 
precise experimental data are now available. For these reasons it is necessary the 
reanalysis of QCD at low energies with account of higher order terms of PT and 
OPE in comparison with experimental data. In this report I present the results of 
such analysis. 



2 Test of QCD at low energies on the basis of 
r-decay data. Determination of ^^(m^) and of 
condensate values. 

Recently, collaborations ALEPH [4], OPAL [5] and CLEO [6] had measured with a 
good accuracy the relative probability of hadronic decays of r-lepton = B{t —>■ 
Vr + hadrons) / B{t z/T-el^e), the vector V and axial A spectral functions and the 
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cross section of e+e~ — > hadrons. Below I present the results of the theoretical 
analysis of r-dccay data basing on the operator product expansion (OPE) in QCD 
[7,8]. In the perturbation theory series the terms up to ctg will be taken into account, 
in OPE - the operators up to dimension 8. 

Consider the polarization operator of hadronic currents 

- ^ / e^'^%TJ,{x)MQ)^)dx = {q,q. - q,.q')Il'-j\q') + q,qM^\q'), (1) 

where J ^V,A; V^^ uj^d, = u'y^-f^d. 
The spectral functions measured in r-decay are imaginary parts of nj^(s) and 

vi/ai{s)= 27r ImU'yj^ (s + iO), oq (s) = 27r ImUf (s + iO) (2) 

Functions Ily\q'^) and 11^^ (g^) are analytical functions in the complex plane 
with a cut along the right-hand semiaxis starting from 4m^ for Ily\q'^) and 9m^ 
for Il^^\q'^). Function 11^^ (g^) has kinematical pole at = 0. This is a specific 
feature of QCD following from chiral symmetry within massless u and d quarks and 
from its spontaneous violation. The kinematical pole appears due to one-pion state 
contribution into n^(g), which has the form [7] 

^ q^ q'^ q^ — 

Consider first the ratio of the total probabihty of hadronic decays of r-lepons into 
states with zero strangeness to the probability of r — > i/T-eVg. This ratio is given by 
the equahty [9] 

B{t + hadronss=o) 

<-: 

(4) 



R. 



T,V+A — 



ud\ i^EW 








1 + 2^ I {vi + ai + ao)(s) - 2^ao(s 



) ml 



where \Vud\ = 0.9735 ± 0.0008 is the matrix element of the Kabayashi-Maskawa 
matrix. Sew = 1.0194 ± 0.0040 is the electroweak correction [10]. Only one-pion 
state is practically contributing to the last term in (4) and it appears to be small: 

Ai?(o) = -24:71^1^ = -0.008 (5) 

Denote 

u{s) =-^1 + 01 + 00 = 2nlm[n^y'^ (s) + n(l)^(s) + (s)] = 2TTlmU{s) (6) 

As follows from eq.(3), n(s) has no kinematical pole, but only right-hand cut. It is 
convenient to transform the integral in eq.(4) into that over the circle of radius 
in the complex s plane [11]-[13]: 
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Calculate first the perturbative contribution into eq.(7). To this end, use the Adler 
function D{Q^): 

D{Q') ^ -2^ '^^^ = E ^na" , a ^ ^ , ^ -s, (8) 

^ ^ n>0 



the perturbative expansion of which is known up to terms ~ Q;^- In MS regulariza- 
tion scheme Ko = Ki = 1, K2 = 1.64 [14], = 6.37 [15] for 3 flavours and for K4 
there is the estimate K4 = 25 ± 25 [16]. The renormgroup equation yields 

da ^, , ^ ^ „,n , r<Q^) da , , 



in the MS scheme for three flavours /3o = 9/4, /3i = 4, ^2 = 10.06, /^a = 47.23 
[17,18]. Integrating over eq.(8) and using eq.(19) we get 

Put /x^ = and choose some (arbitrary) value a{m'^). With the help of eq.(9) 
one may determine then a(<5^) for any and by analytical continuation for any s 
in the complex plane. Then, calculating (10) flnd n(s) in the whole complex plane. 
Substitution of n(s) into eq.(7) determines Rr for the given a(m^) up to power 
corrections. Thereby, knowing Rr from experiment it is possible to flnd the corre- 
sponding to it a(m^). Note, that with such an approach there is no need to expand 
the nominator in eqs.(9),(10) in the inverse powers of InQ'^/fi'^. Particularly, there 
is no expansion on the right-hand semiaxis in powers of the parameter tt / ln{Q'^ / /j,'^) , 
which is not small in the investigated region of Q^. Advantages of transformation of 
the integral over the real axis (4) in the contour integral are the following. It can be 
expected that the applicability region of the theory presented as perturbation the- 
ory (PT) -|- operator expansion (OPE) in the complex s-plane is off the shadowed 
region in Fig.l. It is evident that at positive and comparatively small s PT-I-OPE 
do not work. At negative s — —Q^ in as order a nonphysical pole appears, in higher 
orders, according with (9) it is replaced by a nonphysical cut, which starts from the 
point —Ql, determined by the formula 

ln^ = - r (11) 
//^ Ja{n'^) /3{a) 

Integration over the contour allows one to obviate the dashed region in Fig.l (except 
for the vicinity of the positive semiaxis, the contribution of which, is suppressed by 
the factor (1 — ^)^ in eq.(7)), i.e. to work in the applicability region of PT-I-OPE. 
The OPE terms, i.e., power corrections to polarization operator, are given by the 
formula: 
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Figure 1: The applicability region of PT and OE in the complex plane s. In the 
dashed region PT + OE does not work. 
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(a^-corrections to the 1-st and 2-d terms in eq.(ll) were calculated in [19] and [20], 
respectively). Contributions of the operator with d = 2 proportional to m^, 
and of the condensate 2(m„ + md){0\qq\0) are neglected. (The latter is of an order 
of magnitude smaller than the gluonic condensate contribution) . When calculating 
the d=6 term, factorization hypothesis was used. It can be readily seen that d=4 
condensates (up to small ag corrections) give no contribution into the integral over 
contour cq.(7). The contribution from the condensate (Og) may be estimated as 
I (Og) I < 10~^GeV^ and appears to be negligibly small. Rry+A may be represented 
as 

Rr,v+A = 3|Kd|'^w (l + S'ew + S^''^ + S'^l^) + Ai?(°) = 

= 3.475 ±0.022 (13) 

where = (5/127r)Q;em("^r) = 0.001 is electromagnetic correction [21], Sj^^y — 
— (5 ± 2) • 10^^ is the contribution of d=6 condensate (see below) and 6^^^ is the 
PT correction. The right-hand part presents the experimental value obtained as 
a difference between the total probability of hadronic decays — 3.636 ± 0.021 
[22] and the probability of decays in states with the strangeness S — —1 Rr,s — 
0.161 ± 0.007 [23,24]. For perturbative correction it follows from eq.(13) 

5^°) = 0.206 ±0.010 (14) 

Employing the above described method in ref.[8] the constant Q;s(m^) was found 
from (14) 

a,(m2) = 0.355 ±0.025 (15) 

The calculation was made with the account of terms ~ a^, the estimate of the effect 
of the terms ~ is accounted for in the error. May be, the error is underesti- 
mated (by 0.010-0.015), since the theoretical and experimental errors were added in 
quadratures. 

I determine now the values of condensates basing on the data [3]- [4] on spectral 
functions. It is convenient first to consider the difference Ily — 11^, which is not 
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contributed by perturbative terms and there remains only the OPE contribution: 



- n?(^) = E ( 1 + ) (16) 



D>4 



0^/ 



The gluonic condensates contribution falls out in the V — A difference and only the 
following condensates with d=4,6,8 remain 

O^-^ ^2{mu + md) <qq>= - />' (17) 
O^-^ = 27ra, ( (iZ7^A«d)(J7^A«ii) - {u^5l^.X''d)id^,^^X''u) ) = 

= - — ^ < m >' (18) 

A" 

Og^-^ = 87ra, < gg >2 , - g(0 | qa^^ — G;,q \0)=ml{0\qq\ 0) (19) 

In the right-hand of (18) and the first of eq.'s (19) the factorization hypothesis 
was used. Calculation of the coefficients at ag in eq.(16) gave C4 = 4/3 [19] and 
C6 = 89/48 [20]. The value of ^^(m^) (15) corresponds to as{lGeV^) = 0.60. 
Thus, if we take for quark condensate at the normlization point /J,^ — IGeV^ the 
value following from Gcll-Mann-Oakes-Renner relation at m„ + rrid = 11.7 MeV, 
then vacuum condensates with the account of cts-corrections appear to be equal (at 
//^ = IGeV^y. 

O4 = -4.22 • 10-^ GeV"^ (20) 

Oe = -3.75 • 10-^ GeV^ (21) 

Os = 2.5 ■ 10--'^ GeV^ (22) 

(In what follows, indeces V — A will be omitted and Od will mean condensates with 
the account of corrections) . 

Our aim is to compare OPE theoretical predictions with experimental data on 
V — A structure functions measured in r- decay and the values of Oq and Og found 
from experiment to compare with eqs.(21),(22). Numerical values of Og and Og 
(21), (22) do not strongly differ. This indicates that OPE asymptotic series (16) at 

= — s ~ IGeV^ converge badly and, may be, even diverge and the role of higher 
dimension operators may be essential. Therefore it is necessary to improve the series 
convergence. The most plausible method is to use Borel transformation. Write for 
Hy^ — n^'' the subtractionless dispersion relation 

Put s — s^Q {(f) — on the upper edge of the cut) and make the Borel transformation 
in So- As a result, we get the following sum rules for the real and imaginary parts 
of (23): 



2fc+2 



(24) 



5 




Figure 2: Eq.(24): the left-hand part is obtained basing on the experimental data, 
the shaded region corresponds to experimental errors; the right-hand part - the 
theoretical one - is represented by the dotted curve, numerical values of condensates 
are taken to be equal to the central values of eqs.(26),(27); a) ^ = 37r/4, h) (f) — 57r/6. 



r f ■'^ C\ ■ f ''^ ■ A\f \f ds ^ fe Sin(fc0) Q2fc+2 

(25) 

The use of the Borel transformation along the rays in the complex plane has a 
number of advantages. The exponent index is negative at 7r/2 < < 37r/2. Choose 
in the region 7r/2 < < tt. In this region, on one hand, the shadowed area in 
fig.l in the integrals (24), (25) is touched to a less degree, and on the other hand, 
the contribution of large s, particularly, s > , where experimental data are 
absent, is exponentially suppressed. At definite values of the contribution of some 
condensates vanishes, what may be also used. In particular, the condensate Og does 
not contribute to (24) at = 57r/6 and to (25) at = 27r/3, while the contribution 
of Oe to (24) vanishes at = 37r/4. Finally, a well known advantage of the Borel sum 
rules is factorial suppression of higher dimension terms of OPE. Figs. 2, 3 presents 
the results of the calculations of left-hand parts of eqs.(24),(25) on the basis of the 
ALEPH [1] experimental data comparing with OPE predictions - the right-hand 
part of these equations. 

The experimental data are best described at the values [7] 

Oe = -(6.8±2.1) • 10-=^ Gey^ (26) 

Og = -(7 ± 4) • 10-^ GeV^ (27) 

When estimating errors in (26), (27), an uncertainty of higher dimension operator 
contribution was taken into account in addition to experimental errors. (For details 
-see [7]). 

As is seen from the figures, at these values of condensates a good agreemeent 
with experiment starts rather early - at > Q.bGeV^. The values (26), (27) are 
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Figure 3: The same for eq.(53): a) = 27r/3, b) = 37r/4. 

by a factor of 1.5-2 larger than (21), (22). As was discussed above, the accuracy of 
(21), (22) is of order 50%. Therefore, the most plausible is that the real value of 
condensates Oq, Og is somewhere close to the lower edge of errors in (26), (27). 

Consider now the polarization operator n(s) defined in (6) and condensates en- 
tering OPE for n(s) (see (12)). In principle, the perturbativc terms contribute to 
chirality conserving condensates. If we will follow the separation method of pertur- 
bative and nonperturbative contribution by introducing infrared cut-off [25,26], then 
such a contribution would really appear due to the region of virtualities smaller than 
11^ . In the present paper, according to [8], an another method is exploited, when the 
/9-function is expanded only in the number of loops, (see eq.(lO) and the text after 
it) but not in 1/lnQ^. So, the dependence of condensates on the normalization point 
//^ is determined only by perturbative corrections, as is seen in (12). Condensates 
determined in such a way may be called n-loop ones (in the given case - 3-loop). 
Consider the Borel transformation of the sum Il{s)pert + ^{s)nonpert where Il{s)pert 
is given by eq.(lO), and Il{s)nonpert - by eq.(12). Fig. 4 presents the results of 3-loop 
calculation for two values of Q;s(m^) - 0.355 and 0.330. The widths of the bands 
correspond to theoretical error taken to be equal to the last accounted term K^a^ 
in the Adler function (8). (The same result for the error is obtained if one takes 4 
loops in /9-function and puts K4 = 50 ±50). The dotted line corresponds to the sum 
of contributions of gluonic condensate (^G^) = 0.012 GeV^ and O^^^ condensate 
in (12) with numerical value corresponding to Oq~^ (26). The dots with errors 
present experimental data. (The contribution of the operators d — A and d — 6 is 
given separately in the insert). 

It is seen that the curve with tt,s(m^) = 0.330 and condensate contributions 
can be agreed with experiment, starting from = lAGeV"^, the agreement being 
improved at smaller values (0|^G^|0) than 0.012 GeV*^. The curve with ^^(m^) = 

0. 355 with the account of condensates coincides with experiment only at > 

1. bGeV'^. The same tendency persist for the Borel sum rules taken along the rays 
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M^, GeV2; 

.l I .... I . I. .. I .... I .... I .... I .... I .... I . I. .. I ... . 
0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 

Figure 4: The results of the Borel transformation of V + A correlator for two values 
asiiTT-T = 0.355 and ^^(m^) = 0.330. The widths of the bands correspond to PT 
errors, dots with errors - experimental data. The dotted curve is the sum of the 
perturbative contribution at as{ml) = 0.330 and O4, Oq condensates. 

in the s complex plane at various (p. Fig.5 gives the sum rule for (p — Stt/G. Prom 
consideration of this and of other sum rules there follows the estimation for gluonic 
condensate: 

(0 I -^Gl^G^^ I 0) = 0.006 ± 0.012 GeV^ (28) 

The best agreement of the theory with experiment in the low region (up to 
~ 2% at Af^ > 0.8Gel^^) is obtained at Q;s(m^) = 0.330 which corresponds to 
a,(m2) = 0.118. 

It was shown [8], that in the dilute instanton gas appoximation [27] instantons 
do not practically affect determination of Q;s(m^) and the Borel sum rules. Their 

effect, however, appears to be considerable and strongly dependent on the value of 
the instanton radius in the sum rules obtained by integration over closed contours 
in the complex plane s at the radii of the contours s < 2GeV^. 

3 Sum rules for charmonium and gluonic conden- 
sate. 

In this Section the charmonium sum rules are revisited. (In what follows I formulate 
the main results of [28]). 

Consider the polarization operator of charmed vector currents 

1 1 dx e*^" ( TJ^{x)MO) ) = ( q^q, - g^^q^ ) U{q'') , = c^^c (29) 
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Figure 5: The result of the Borel transformation along the ray at = Svr/G. The 
dotted line corresponds to the central value of gluonic condensate (56), added to the 
PT curve with as{ml) = 0.330. 



The dispersion representation for n(5^) has the form 

R{s) = 47rlmn(. + iO) , n(g^) = H , (30) 

where R{oo) = 1 in partonic model. In approximation of infinitely narrow widths 
of resonances R{s) can be written as sums of contributions from resonances and 
continuum 

3 71 

^(*) = n-2 2 ( \ Y.rn^^^^eeS{s - ml) + 9{s - Sq) (31) 

where Qc = 2/3 is the charge of charmed quarks, sq - is the continuum threshold 
(in what follows ^/sq = 4.6GeV), a{s) - is the running electromagnetic con- 
stant, (^{iTi^^^) = 1/133.6 Following [1], to suppress the contribution of higher 
states and continuum we will study the polarization operator moments 

According to (31) the experimental values of moments are determind by the equality 

It is reasonable to consider the ratios of moments M„i((5^)/M„2(Q^) from which 
the uncertainty due to error in Fj/^^ee markedly falls out. Theoretical value for 



9 



— I — I — I — I — I — I — I — I — I — I — 1— 



1.32 



1.3 



1.28 



1.26 



1.24 



m, GeV 



— I — I — I — I — I — I — I — I — I — I — r- 




-0.03 -0.02 -0.01 0.01 0.02 0.03 



Figure 6: The dependence of m{m) on (0|as/7r)G'^|0) obtained eA n — 10, Q'^ 
0.98 • and a^iQ^ + m^). 



n(g^) is represented as a sum of perturbative and nonperturbative contributions. It 
is convenient to express the perturbative contribution through R{s), making use of 
(30), (32): 

Ris) = ^i?(")(5,/^')«"(/^') (34) 

n>0 

where a(/i^) = as{^'^)/T^- Nowadays, three terms of expansion in (34) are known: 
[29] R^^^ [30], i?*^^) [31]. They are represented as functions of quark velocity 



V = \Jl — Aw? / s, where m - is the pole mass of quark. Since they are cumbersome, 
I will not present them here. 

Nonperturbative contributions into polarization operator have the form (re- 
stricted by d=6 operators): 

^n^,ert{Q') = (^(0 I I 0) [ f^'\z) + af^'\z) ] + 

+ J^^''^^'^^^ I I 0)F(.), . = (35) 

Functions f'^'^\z), f'^^z) and F{z) were calculated in [1], [32], [33], respectively. 
The use of the quark pole mass is, however, inacceptable. So, it is reasonable to 
turn to MS mass m(/x^), taken at the point jj? — fn?. After turning to the MS 
mass m(m^) we 

m(2) m(^'1) 
n=10,Q^ = 4m^: jj^ = 0M5. jp^ = 1-136 , ^ = -1.673 

(36) 

At a ~ 0.1 and at the ratios of moments given by (36) there is a good reason to 
believe that the FT series well converges. Such a good convergence holds (at n > 5) 
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Figure 7: The dependence of m{m) on (0|(as/7r)G'^|0) obtained from the moments 
(horizontal bands) and their ratios (vertical bands) at different a^. The left-hand 
figure: Q"^ = Arff, n = 10, M10/M12; the right-hand figure - = 8m'^, n = 15, 
M15/M17. 

only in the case of large enough Q'^, ai Q'^ — one does not succeed in finding such 
n, that perturbative corrections, corrections to gluonic condensates and the term 
~ (G^) contribution would be simultaneously small. 

It is also necessary to choose the scale - normalization point /j,^ where Q;s(/i^) is 
taken. In (34) R{s) is a physical value and cannot depend on /i^. Since, however, 
we take into account in (34) only three terms, at unsuitable choice of //^ such /i^ 
dependence may arise due to neglected terms. At large the natural choice is 

= Q^. It can be thought that at = the reasonable scale is /i^ = m^, though 
some numerical factor is not excluded in this equality. That is why it is reasonable 
to take interpolation form 

A^' = g' + m^ (37) 

but to check the dependence of final results on a possible factor at frf. Equalling 
theoretical value of some moment at fixed Q'^ (in the region where M^'^ and M^^^ 
are small) to its experimental value one can find the dependence of m on ((a^/Tr)^^) 
(neglecting the terms ~ {G^))- Such a dependence for n ~ 10 and Q"^ /im^ — 0.98 
is presented in Fig. 6. 

To fix both m and {{as/'n)G'^) one should, except for moments, take their ratios. 
Fig.7 shows the value of m obtained from the moment Mio and the ratio M10/M12 
at — Am? and from the moment Afis and the ratio M15/M17 at = Sm?. The 
best values of masses of charmed quark and gluonic condensate are obtained from 
Fig.7: 

m(m2) = 1.275 ±0.015 GeV , ^^G^^ = 0.009 ± 0.007 GeV^ (38) 
Up to now the corrections ~ (G^) were not taken into account. It appears that in 
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the region of n and used to find m and gluonic condensate they are comparatively 
small and, practically, not changing rn. increase {{as/7i)G'^) by 10 — 20% if the term 
~ (G^) is estimated according to instanton gas model [34] at pc = 0.5/m. 

It should be noted that improvement of the accuracy of Fj/^^ee would make it 
possible to precise the value of gluonic condensate: the widths of horizontal bands 
in fig.7 are determined mainly just by this error. In particular, this, perhaps, would 
allow one to exclude the zero value of gluonic condensate, that would be extremely 
important. Unfortunately, eq.(38) does not allow one to do it for sure. Diminution 
of theoretical error which determine the width of vertical bands seems to be less 
real. 

4 Conclusion 

In this report I compare the results of the recent precise measurements of r-lepton 
hadronic decays [4]-[6] with QCD predictions in the low energy region. The pertur- 
bative terms up to and the terms of the operator product expansion (OPE) up to 
d=8 were taken into account. It is shown that QCD with the account of OPE terms 
agrees with experiment up to ~ 2% at the values of the complex Borel parameter 
|M^| > 0.8 — l.OGeV^^ in the left-hand semiplane of the complex plane. It was found: 

1. The values of the QCD coupling constant a;s(m^) = 0.355 ± 0.025 from the total 
probability of r-decays and a.,(m^) = 0.330 from the sum rules at low energies. 
(The latter value corresponds to ^^(m^) = 0.118). 

2. The value of the quark condensate square (assuming factorization) 

a,{\ i^i/j I 0)2 = (2.25 ± 0.70) • 10"^ GeV^ 

and of quark-gluon condensate of d=8. 

3. The value of gluonic condensate: 

a) from the r-decay data: 

(0 I —G^ I 0) = (0.006 ± 0.012) GeV^ 

TT 

b) from the sum rules for charmonium 

(0 I —G^ I 0) = (0.009 ± 0.007) GeV^ 

71 

It is shown that the sum rules for charmonium are in agreement with experiment 
when accounting for perturbative corrections ~ and for OPE terms proportional 
to {{as/7r)G^) and to (G^). 

The main conclusion is that in the range of low-energy phenomena under con- 
sideration, perturbation theory and operator product expansion are in an excellent 
agreement with experiment starting from ~ IGeV^. 

I am deeply indebted to K.N.Zyablyuk who had made the main calculations in 
papers [7,8,28], the results of which I used here. 

This work was supported by the grants CRDF RP2-2247, INTAS-2000-587 and 
RFFI 00-02-17808. 
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